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Recently, it was shown that differential rotation is an unavoidable feature of nonlinear r-modes. 
We investigate the influence of this differential rotation on the detectability of gravitational waves 
emitted by a newly born, hot, rapidly-rotating neutron star, as it spins down due to the r-mode 
instability. We conclude that gravitational radiation may be detected by the advanced laser inter- 
ferometer detector LIGO if the amount of differential rotation at the time the r-mode instability 
becomes active is not very high. 
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I. INTRODUCTION 

The new kilometer-scale laser interferometer 
gravitational-wave detectors LIGO and Virgo, or 
their advanced versions, operating at a broad frequency 
band (between about 10 and 10^ Hz), may detect 
gravitational waves from a variety of sources (for a 
recent review see Ref. Q). One such source could be a 
spinning neutron star emitting gravitational radiation 
due to the r-mode instability. 

First studied more than twenty years ago , r-modes 
are non-radial oscillation modes of rotating stars that 
have the Coriolis force as their restoring force and a char- 
acteristic frequency comparable to the angular velocity 
of the star. These modes have been attracting increasing 
attention since it was discovered that they are driven un- 
stable by gravitational radiation emission 0,Q and that, 
for a large range of relevant temperatures and angular 
velocities of neutron stars, the driving effect of gravita- 
tional radiation is stronger than the damping effect of 
viscosity 

Using a phenomenological model for the evolution of 
the r-mode instability, Owen et al. 7J studied the de- 
tectability of gravitational waves emitted by newly born, 
hot, rapidly-rotating neutron stars, arriving at the con- 
clusion that such waves could be detected by the en- 
hanced version of LIGO if sources were located at dis- 
tances up to 20 Mpc from Earth. However, a deeper 
understanding of this issue requires taking into account 
nonlinear effects in the evolution of the r-mode instabil- 
ity. Arras et al. Q considered the nonlinear interaction 
between modes, arriving at the conclusion that enhanced 
laser interferometer detectors could see gravitational ra- 
diation from r-modes if the sources were much closer to 
Earth, namely, at distances smaller than about 200 kpc. 
In this paper we will investigate the issue of detectability 
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of gravitational radiation emitted by newly born neutron 
stars as they spin down due to the r-mode instability, but 
taking into account another nonlinear effect: differential 
rotation. 

Rezzolla et al. were the first to suggest that r-modes 
induce a drift of fluid elements along azimuthal direc- 
tions and derived an approximate analytical expression 
for these drifts. Soon afterwards, numerical studies, both 
in gene ral relativistic [lfil | and Newtonian hydrodynam- 
ics jij, confirmed the existence of such drifts. Recently, 
an exact r-mode solution, representing differential rota- 
tion of kinematic nature that produce large scale drifts 
of fluid elements along stellar latitudes, was found within 
the nonlinear theory up to second order in the mode's 
amplitude in the case of a Newtonian, barotropic, non- 
magnetized, perfect- fluid star This differential ro- 
tation, which is an unavoidable feature of nonlinear r- 
modes, contributes to the physical angular momentum 
of the mode [l^ and, therefore, plays an important role 
in the nonlinear evolution of the r-mode instability ^3 ■ 
In particular, the amplitude of the r-mode saturates in 
a natural way at a value that depends on the amount of 
differential rotation at the time the instability becomes 
active; this saturation value can be much smaller than 
unity 13]. 

In this paper, using the model of evolution of Ref. [l3l |. 
we investigate how differential rotation induced by r- 
modes influences the detectability, by the laser interfer- 
ometer detectors LIGO and Virgo, of gravitational waves 
emitted by a newly born, hot, rapidly-rotating neutron 
star. In Sect. ^ we briefly review the evolution model 
of the r-mode instability. The main results of the paper 
are presented in Sect. IIIII where we derive expressions 
for the maximum value of the gravitational wave ampli- 
tude and for the frequency-domain gravitational wave- 
form. We also compare, for different values of differential 
rotation, the characteristic amplitude of the signal with 
the rms strain noise in the initial LIGO, Virgo and ad- 
vanced LIGO gravitational-wave detectors and compute 
the signal-to-noise ratio. Finally, Sect. lIVl is devoted to 
a discussion of the results and conclusions. 
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II. NONLINEAR EVOLUTION OF THE 
ii-MODE INSTABILITY 

In a recent paper 0, we have investigated the role 
of differential rotation in the evolution of the I = 2 r- 
mode instability in a newly born, hot, rapidly- rotating 
neutron star, using a simple phenonienological model 
adapted from the one proposed in Rcf. 7J. The main 
difference between our modified model and the one pro- 
posed in Ref. 't'I is that the former takes into account that 
the full physical angular momentum of the r-mode per- 
turbation also includes a contribution from differential 
rotation. Namely, at second order in the mode's ampli- 
tude a, the physical angular momentum was taken to be 

El 



(1) 



where il, R and M denote the angular velocity, the ra- 
dius and the mass of the star, respectively, and J = 
pr^dr/{MR*) = 1.635 x IQ-^. In the above expres- 
sion, it was assumed that the star's mass density p and 
pressure p are related by a polytropic equation of state 
p = kp^ with k such that M = IAMq and R = 12.53 km. 
In Eq. (Q, K IS a constant fixed by initial data, giving 
the initial amount of differential rotation associated with 
the r-mode. 

Within our model (see Ref. for details), the total 
angular momentum is given by the sum of the angular 
momentum of the unperturbed star and the angular mo- 
mentum of the r-mode perturbation. 



(2) 



where / = (Stt/S) Jq pr^dr = IMR^ (/ = 0.261) is the 
momentum of inertia of the unperturbed star. Assuming 
that the total angular momentum of the star decreases 
due to the emission of gravitational radiation and that 
the angular momentum of the perturbation increases due 
to the emission of gravitational radiation and decreases 
due to the dissipative effect of viscosity, we arrive at a 
system of differential equations determining the evolu- 
tion of the star's angular velocity ri(<) and the r-mode's 
amplitude a(t), namely. 
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where Q = 3 J/(2/) = 0.094 and the value of K is chosen 
to lie in the interval —5/4 ^ A' ^ 10^^. The upper limit 
for K results from the fact that one wishes to impose the 
condition that the initial value of the angular momentum 
of the r-mode is much smaller than the angular momen- 
tum of the unperturbed star, i.e., S^^^Jq <C I^o implies 



that K <€. 10^^ if we choose ao = 10~^. The lower Hmit 
for K results from the fact that we do not want to sat- 
urate the amplitude of the mode by hand, a procedure 
needed for the case K < —5/4 in order to avoid that the 
total angular momentum of the star becomes negative 
(see Ref. 0| for details). 

In Eqs. © and 10} the gravitational-radiation and vis- 
cous timescales^ are given, respectively, by ^ 
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with the fiducial timescales tgr, ~ —3.26 s, — 2.52 x 
10^ s and TB = 2.01 x 10" s ,15|. The temperature 
of the star is assumed to decrease due to the emission of 
neutrinos via a modified URCA process. 
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where Tq is the initial temperature of the star and = 
1 yr characterizes the cooling rate • 

The viscous timescale given above was derived for a 
simple model of a neutron star with shear and bulk vis- 
cosity. The consideration of other types of viscosity (as, 
for instance, viscosity at the core-crust boundary in neu- 
tron stars with a crust or hyperon bulk viscosity TT]) 
leads to expressions for the viscous timescale that differ 
substantially from Eq. However, because of the com- 
plexity of the involved processes it is not clear yet which 
viscous timescale is more appropriate to describe a real 
neutron star. In this paper we choose to restrict ourselves 
to the viscous timescale given by Eq. ©. 

Let us assume that the temperature, the mode's am- 
plitude and the star's angular velocity take the initial 
values To = 10"'^^ K, ao = 10~^ and fig = ^k, where 
ilif — {2/3)^/TrGp = 5612 s^^ is the Keplerian angu- 
lar velocity at which the star starts shedding mass at 
the equator (these initial values will be used throughout 
the paper). Then, in the first moments of the evolu- 
tion, bulk viscosity dominates the dynamics of Eqs. 
and 10}. But this lasts just a fraction of a second. In- 
deed, the temperature of the star falls so rapidly that the 
gravitational-radiation driving effect almost immediately 
dominates over the bulk viscosity damping effect. This 



^ The gravitational-radiation timescale tqh for r-modes was first 
obtained using an expression for the time evolution of the physi- 
cal energy of the r-mode perturbation, assuming that the imag- 
inary part of the frequency of the mode, Imiuj) = I/tqh, is re- 
lated to the time derivative of the energy by dE/dt = —2E/tqii 
13 . But TQji can also be obtained by solving explicitly the hydro- 
dynamic equat ion s in the presence of the gravitational radiation 
reaction force 
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preponderance of gravitational radiation continues dur- 
ing most of the evolution, ending between = 3.6 x 10^ s 
(for K ^ -5/4) and U ^ 7.1 x 10^ s (for iiT > 1) . Af- 
terwards, shear viscosity dominates the dynamics of the 
evolution. Ki t — 1 yr, which corresponds to a temper- 
ature T = 10^ K, we stop evaluating Eqs. © and Q), 
since at this temperature superfluid effects are expected 
to become important, rendering invalid our assumptions 
about the viscous timescale [l8l |. 

In Fig. ^ the numerical solution of the system of equa- 
tions ||2J) and (01, for different values of the constant is 
represented in a (fi, T) diagram. As we can see there, for 
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As shown in Ref. Ts'l, Eqs. © and © yield the fol- 
lowing solution: 
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An analysis of the above solution makes it clear that 
the evolution of the r-mode instability proceeds in two 
stages. 

In the first stage of the evolution, the above solution 
is well approximated by 

^{i) , 4,^^ „ { f \^ 2t 



FIG. 1: Angular velocity of the star as a function of the 
temperature. As the star cools and emits gravitational radi- 
ation, its angular velocity decreases to values which are quite 
insensitive to the value of K. The dotted line represents the 
stability curve, TQlf^iyt) -\-Ty^{yL, T) = 0, i.e., the set of points 
for which the damping effect of viscosity balances exactly the 
driving effect of gravitational radiation. 

a newly born, hot, rapidly-rotating neutron star, there is 
an interval of relevant temperatures and angular veloci- 
ties of the star for which the r-mode instability is active. 
For such values of temperature and angular velocity (cor- 
responding to < i < ff,), the influence of viscosity on 
the evolutionary equations (O and l@J is so small, as 
compared with the influence of gravitational radiation, 
that n,{t) and a[t) can be determined, in a very good 
approximation, by the system of equations: 
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Therefore, in what follows, Eqs. ^ and lO, which can 
be solved analytically, are used to determine fl(t) and 
a{t), while Eqs. Q and which we solve numerically, 
are used just to determine tb, i.e., the time at which the 
damping effect of viscosity becomes equal to the driving 
effect of gravitational radiation. As mentioned above, 
the value of tf,, which depends on lies in the range 
3.6 X 10^ s 4 7.1 X 10^ s. 
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i.e., both the angular velocity of the star and the am- 
plitude of the mode evolve on the gravitational-radiation 
timescale. 

In the second stage of evolution, the solution (|10|) and 
pi|l is well approximated by 
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i.e., due to nonlinear effects (the presence of differen- 
tial rotation) , the amplitude of the mode saturates at a 
value that depends crucially on the parameter K, namely, 
Usat oi {K + 2)"-^/^, and the angular velocity decreases 
to a final value which is quite insensitive to the value of 
K (it depends on K through tf,). 

The smooth transition between the two stages of evo- 
lution occurs a few hundred seconds after the instabil- 
ity becomes active. The moment of occurrence of this 
transition can be defined more precisely by introducing 
the transition time ta, corresponding to the moment, de- 
termined by the condition d'^aj dt^{ta) — 0, when the 
mode's amplitude changes from an exponential to a much 
slower power-law growth. Using Eqs. (jSJ, and (O, the 
above condition yields a{ta) = [\2{K+2)Q]^'^/'^ or, using 
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the relation 17 « no[l + ^{K + 2) Qa'^]-\ n{ta) = 0.9f7o. 
Inserting a{ta) into Eq. imj, or fl{ta) into Eq. ((Tn)l . 
and taking into account that 4:{K + 2)QaQ/3 <C 1 and 
ao ^ a(ia), one obtains 



ta « [521.0- 18.5 ln(ii: + 2)] s. 



(16) 



III. DETECTABILITY OF GRAVITATIONAL 
WAVES FROM THE fl-MODE INSTABILITY 

The gravitational radiation emitted by a newly born 
neutron star, as it spins down due to the r-mode insta- 
bility, could in principle be detected by the laser inter- 
ferometer detectors LIGO and Virgo using appropriate 
detecting strategies. 

The gravitational wave amplitude, averaged over 
source and detector orientation, is given by 



m\ = 1.3 X 10-"a(f) 



D 



where D is the distance to the source. 

In Fig. 121 the gravitational wave amplitude is shown 
for different values of K and a distance to the source 
of D = 20 Mpc. During the first stage of evolu- 
tion, \h{t)\ grows exponentially [see Eqs. ((T^ and ifT^ ]. 
while during the second stage of evolution it decreases 
slowly, as i'^/^ [see Eqs. (O and US))]. The maximum 
value h„iax = |^(^*)| is achieved when d\h\/dt{t^) = 0, 
which, using Eqs. © and ||5J), gives the mode's ampli- 
tude Q!(i*) = \/3/[20{K + 2)Q], or, using the relation 
ft w ilo[l + ^{K + 2)Qa'^]~^, the star's angular velocity 
ri(i*) = {5/6)ilK- Substituting these values of a and H. 
into Eq. l(T7|l . we obtain 
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If Q!(t*) is inserted into Eq. lfTT|l . or f7(i*) into Eq. ifTUIl .we 
obtain that the gravitational wave amplitude reaches its 
maximum value at ~ [543.1 — 18.51n(if -f 2)] s, i.e., 
shortly after the transition between the first and second 
stages of evolution of the r-mode instability. 

Let us now derive an analytical expression for the 
frequency-domain gravitational wave amplitude [/i(/)[. 
Using the stationary phase approximation, namely, 
\h{f)\ = \h{t)\/^\df/dtl where df/dt is straightfor- 
wardly obtained from Eq. ||SJ|, 
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we can easily obtain the frequency-domain gravitational 
wave amplitude \h{f)\, 
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FIG. 2: Gravitational wave amplitude \h{t)\ for different 
values of K and distance to the source D = 20 Mpc. At 
t, ~ [543.1 — 18.5 ln(A'+2)] s the gravitational wave amplitude 
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where / — 2il/(37r) is the frequency of the emitted grav- 
itational wave. 

Equation (|20|l applies to both stages of evolution, 
for frequencies ranging from fmin — (77 — 80) Hz to 
fmax = 1191 Hz, where the maximum frequency cor- 
responds to the initial value of the angular velocity of 
the star, fix = 5612 s""'^, and the minimum frequency 
corresponds to the final angular velocity of the star, 
il(tb), which lies between 0.06517;^ (for K = —5/4) and 
OMTQk (for if > 1). 

Let us mention that in Ref. |3| an expression for the 
frequency-domain gravitational wave amplitude [/i(/)[ 
was derived based on the assumption that dJ/df oc /, 
where J is the total angular momentum of the star, / 
is the frequency of the emitted gravitational wave and / 
is the moment of inertia of the unperturbed star. Since 
within the model of Ref. |^ , the condition dJ/df cx / only 
applies during the second stage of evolution, the expres- 
sion obtained there for [/»(/)[ is valid just for this stage 
of evolution. Within our model, however, it can be eas- 
ily shown that the condition dJ/df cx / holds, not only 
during the second stage of evolution, but also during the 
first stage. Indeed, from Eqs. and (0), we obtain that 
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Now we can use Eqs. ((HJ and ^ to relate dil and da, 
namely. 
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Inserting this expression into Eq. H21(l and using / = 
2ri/(37r) we obtain 



dJ 97rJ 



df 8{K + 2)' 



(23) 
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Since the above condition holds during both stages of 
the evolution, we could use it to derive, in the manner 
proposed in Ref. , the frequency-domain gravitational 
wave amplitude given by Eq. (|2U|I . 

The gravitational wave amplitude \h{f)\ is represented 
in Fig. 13 for different values of K and for gravitational- 
wave frequencies lying between fmin and fmax- 
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FIG. 3: Gravitational wave amplitude in the frequency do- 
main, \h{f)\, for different values of K and distance to the 
source D — 20 Mpc. 



It is worth mentioning that within the model of evo- 
lution proposed in Ref. Q, the frequency-domain gravi- 
tational wave amplitude has a spike at high frequencies, 
due to the fact that during the first stage of evolution the 
angular velocity of the star evolves very slowly on the vis- 
cous timescale, leading to a quasi-monochromatic gravi- 
tational wave emission during the first 500 s of evolution. 
However, as we have seen, if one takes into account the 
influence of differential rotation, namely, the fact that it 
contributes to the physical angular momentum of the r- 
modc perturbation, then the angular velocity of the star 
evolves in the gravitational-radiation timescale already 
in the first stage of evolution [see Eq. ©J. As a con- 
sequence, the gravitational wave amplitude |/i(/)| in the 
first stage of evolution is also given by Eq. (|20|1 and we 
observe no spike. 

Let us now analyze the possibility of detecting the 
emitted gravitational waves with laser interferometer de- 
tectors LIGO and Virgo. Because of the complexity of 
real neutron stars, our knowledge of the evolution of the 
r-mode instability is insufficient to predict the gravita- 
tional waveform with such an accuracy that a matched 
filtering signal-processing technique would be feasible. 
However, we can use matched filtering in order to es- 
timate the detectability of the gravitational- wave signal. 
For that purpose, the characteristic amplitude of the sig- 
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is compared with the rms strain noise in the detector, 

hrmsif) = ylKU), (25) 

where Sh{f) is the noise power spectral density of the 
detector. 

For frequencies in the interval 50 Hz ^ / ^ 1200 Hz, 
the curves for the noise power spectral densities of the 
initial LIGO Hi, Virgo 20] and advanced LIGO HI de- 
tectors are well approximated by the following analytical 
expressions, respectively: 



ShU) = s. 



where Si = 3.4 x lO^^^ Hz"^ and /i = 142.0 Hz, 
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where ^2 = 1.5 x lO""*^ Hz"' and /a = 249.6 Hz, and 
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where 5*3 2.2 x lO"''^ Hz'^ /g = 52.8 Hz and f^ = 
421.3 Hz. 

In Fig. ^ the curves corresponding to the rms strain 
noises in the initial LIGO, Virgo and advanced LIGO 
detectors are compared with the curves corresponding 
to the characteristic amplitude of the gravitational wave 
signal for different values of K and for a distance to the 
source of _D = 20 Mpc. For such a distance, which in- 
cludes the Virgo cluster of galaxies, a few supernovae per 
year are expected. 

The most striking feature of Fig. 0| is that the de- 
tectability of gravitational waves from the r-mode insta- 
bility of newly born neutron stars is drastically reduced 
as the initial value of differential rotation associated with 
r-modes increases. Indeed, for K > 100, even the ad- 
vanced LIGO detector would not have enough sensitivity 
to see such sources (for D = 20 Mpc). 

The visual comparison, in Fig. ^ between the charac- 
teristic amplitude of the signal and the rms strain noise in 
the detector gives us a qualitative measure of the signal- 
to-noise ratio for matched filtering. A quantitative de- 
termination of the signal-to-noise ratio is obtained from 
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FIG. 4: The dotted lines correspond to the rms strain noises 
in the initial LIGO, VIRGO and advanced LIGO detectors. 
The solid lines correspond to the characteristic amplitude for 
different values of K. The distance to the source is taken here 
to be D = 20 Mpc. 
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FIG. 5: The dotted line corresponds to the rms strain noise 
in the advanced LIGO detector. The solid lines correspond 
to the characteristic amplitude for different values of K. The 
distance to the source is taken here to be 7? = 30 kpc. 
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For the initial LIGO and Virgo detectors, even for 
small initial differential rotation {K « 0), the signal- 
to-noise ratio is not significant for D = 20 Mpc. Of 
course, this ratio can be increased if we consider sources 
located at smaller distances, but at the cost of decreasing 
the number of expected supernova events per year and, 
hence, the probability of a detection. 

For the advanced LIGO detector the situation im- 
proves: for small values of K, the signal-to-noise ratio 
is significant even for D — 20 Mpc. Therefore, and since 
within such a distance several supernovae per year are 
expected, one could hope that advanced LIGO detectors 
would see gravitational radiation from the r-mode in- 
stability of young neutron stars. However, such hope 
is based on the assumption that neutron stars are born 
with small initial differential rotation associated with r- 
modes. But this may not be the case. And if a neutron 
star is born with substantial differential rotation asso- 
ciated with r-modes (say K > 100), then the emitted 
gravitational waves will not be seen even by advanced 
LIGO. Of course, as already mentioned above, one could 
consider smaller distances to the source, but then the 
number of events per year would decrease, decreasing 
the probability of a detection. In Fig. [3 the character- 
istic amplitude for high values of K is compared with 
the rms strain noise in the advanced LIGO detector for a 
distance to the source of 30 kpc, i.e., within our Galaxy, 
in which about 2 supernovae are expected every hundred 
years f22j . 



The fact that an increase of the initial amount of dif- 
ferential rotation associated with r-modes makes it more 
difficult to detect gravitational waves from these modes 
(for a given distance to the source ) can be easily under- 
stood from angular momentum considerations. Indeed, 
it was shown in Ref. , based on an unpublished gen- 
eral argument of Blandford, that the signal-to-noise ratio 
given by Eq. ^ is well approximated by 



2G AJ 



(31) 



where AJ = Jq — J{tb) is the total amount of angular 
momentum carried away by gravitational waves, J(t) is 
the total angular momentum of the star and Jq ~ li^a is 
the initial angular momentum of the star. Using Eqs. 
and the total amount of angular momentum carried 
away by gravitational waves can be written as 



AJ 



a (tb), (32) 



or, since w f7o[l + l^{K + 2)Qa^]''^, as 



AJ _ 3 
X 4(if + 2) 



(:-^), ,33, 



where n(tb) = (0.065 - O.067)r2o- Thus, for X > 1, only 
a small part of the initial angular momentum of the star 
is carried away by gravitational waves^ (see Fig. EJ and. 



^ As shown in Ref. 13|, for K ^ 1, most of the angular momentum 
of the unperturbed star is transferred to the r-mode perturba- 
tion, due to the fact that the fluid develops a strong differential 
rotation. Indeed, after a few hundred seconds of evolution of 
the r-mode instability, the average differential rotation increases 
rapidly, saturating at high values relatively to the initial angular 
velocity of the star. 
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consequently, detection of these waves becomes a more 
difficult task. 
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FIG. 6: Angular momentum carried away by gravitational 
waves, A J, as a function of time for different values of K. For 
K ^ \, only a small part of the initial angular momentum of 
the star is carried away by gravitational waves. 



IV. CONCLUSIONS 

In this paper we have investigated the influence of dif- 
ferential rotation on the detectability of gravitational ra- 
diation emitted by a newly born, hot, rapidly- rotating 
neutron star, as it spins down due to the r-mode insta- 
bility. 

A model of evolution of the r-mode instability that 
takes into account differential rotation [IS] has been 
used to derive the gravitational wave amplitude \h{t)\ 
and its Fourier transform We have shown that 

the maximum value of the gravitational wave amplitude, 
hmax , depends on the amount of differential rotation at 
the time the r-mode instability becomes active, namely, 
hmax {K -\- 2)~^/^. We have also shown that the 
frequency-domain gravitational wave amplitude, 
has no spike at high frequencies, as opposed to the re- 
sults of Ref. Jj. There, the spike of \h{J)\ is due to 
the fact that during the first stage of evolution of the r- 
mode instability the angular velocity of the star evolves 
very slowly on the viscous timescale, leading to a quasi- 
monochromatic gravitational wave emission during the 
first few hundred seconds of evolution. However, as we 
have seen, if one takes into account the influence of dif- 
ferential rotation, then the angular velocity of the star 



evolves in the gravitational-radiation timescale already 
in the first stage of evolution and, as a consequence, dur- 
ing this stage of evolution the gravitational wave is not 
monochromatic and has no spike. 

We have assumed matched filtering in order to in- 
vestigate the detectability of the gravitational wave sig- 
nal. However, our knowledge of the evolution of the 
r-mode instability is insufficient to predict the gravi- 
tational waveform with such an accuracy that an opti- 
mal matched filtering signal-processing technique would 
be feasible. But non-optimal signal-processing strategies 
could be developed such that they would yield results not 
much different from the ones obtained with matched fil- 
tering (see, for instance, the hierarchical search strategies 
proposed in |24j'). Therefore, our results concerning the 
detectability of gravitational waves from r-modes can be 
considered to be a good approximation. 

Assuming matched filtering, the characteristic ampli- 
tude of the signal hc{f) is compared with the rms strain 
noise in the initial LIGO and Virgo detectors, as well as 
with the expected rms strain noise in the advanced LIGO 
detector [see Figs. 0] and El and Eq. ^U^]. We conclude 
that the detectability of gravitational waves from the r- 
mode instability of newly born neutron stars is drasti- 
cally reduced as the initial value of differential rotation 
associated with r-modes increases. For the initial LIGO 
and Virgo detectors, the signal-to- noise ratio obtained 
with matched filtering for sources located at 20 Mpc is 
smaller than unity even for « 0. For advanced LIGO, 
if neutron stars are born with significant differential ro- 
tation associated with r-modes, then detection of gravi- 
tational waves would be possible only if the distance to 
these neutron stars is considerably smaller than 20 Mpc. 
For instance, UK— 10^, a signal-to- noise ratio greater 
than 10 requires a source located no more than about 80 
kpc away [see Eq. i.e., within a sphere containing 

the Milky Way, the Magellanic Clouds and a few more 
small galaxies, in which just a few supernovae per century 
are expected. However, if the initial value of differential 
rotation associated with r-modes is small {K w 0), then 
S/N > 10 can be obtained for D = 20 Mpc. Since within 
a volume of radius 20 Mpc several supernovae per year 
are expected, it is quite possible that advanced LIGO will 
detect gravitational radiation from the r-mode instability 
of young neutron stars. 
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